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A gray-level image texture descriptors based on fractal dimension estimation is proposed in this work. The 
proposed method estimates the fractal dimension using probability (Voss) method. The descriptors are 
computed applying a multiscale transform to the fractal dimension curves of the texture image. The proposed 
texture descriptor method is evaluated in a classification task of well known benchmark texture datasets. 
The results show the great performance of the proposed method as a tool for texture images analysis and 
characterization. 
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I. INTRODUCTION 

Fractals have played an important role in many areas 
with applications related to computer vision and pat- 
tern recognition^^. This broad use of fractal geometry 
is explained by the flexibility of fractals in representing 
structures usually found in the nature. In such objects, 
we observe a level of details at different scales which are 
described in a straightforward manner by fractals, rather 
than through classical Euclidean geometry. 

Most of fractal-based techniques are derived from frac- 
tal dimension concept. Although this concept was de- 
fined originally only for real real fractal objects, it con- 
tains some properties which turn it into a very interesting 
descriptor for any object of real world. Indeed, the fractal 
dimension measures how the complexity (level of details) 
of an object varies along scales. Such definition corre- 
sponds to an effective and flexible means of quantifying 
the spatial occupation of an object as well as physical 
and visual important aspects which characterize an ob- 
ject, like luminance and roughness, for instance. 

Among these fractal techniques, we may cite 
Mult ifr act als 9 11 , Multiscale Fractal Dimension^™ and 
Fractal Descriptors^^. Here we are focused on the later 
approach which demonstrated the best results in such 
kind of application 21 . The main idea of fractal descrip- 
tors theory is to provide descriptors of an object repre- 
sented in a digital image from the relation between the 
values of fractal dimension taken under different obser- 
vation scales. These values provide a valuable informa- 
tion about the complexity of the object in the sense that 
they capture the degree of details at each scale. In this 
way, fractal descriptors are capable of quantifying impor- 
tant physical characteristics of the structure, like fractal 
dimension, but presenting a sensibly richer information 
than a unique number (fractal dimension). 

Here, we proposed a novel fractal descriptors based 
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on probability fractal dimension. We used the whole 
power-law curve of the dimension and applied a space- 
scale transform to emphasize the multiscale aspect of the 
features. Finally, we test the proposed method over two 
well-known datasets, that is, Brodatz and Vistex, com- 
paring the results with another fractal descriptors ap- 
proach showed irP^I and other conventional texture anal- 
ysis methods. The results demonstrated the power of 
probability descriptors, achieving a more precise classifi- 
cation than other classical techniques. 



II. FRACTAL THEORY 

In the last decades we observe a growing number of 
works applying fractal geometry concepts in the solution 
of a wide range of problemP"^. This interest in frac- 
tals is mainly motivated by the fact that conventional 
Euclidean geometry has severe limitations in providing 
accurate measures of real world objects. 



A. Fractal Dimension 

The first definition of fractal dimension provided is 
the Hausdorff dimension. In this definition, we consider 
a fractal object as being a set of points immersed in a 
topological space. Thus, we can use results from Measure 
Theory to define a measure over this object. This is the 
Hausdorff measure expressed through: 

CO 

H S 5 (X) = inf \ U i\ S '- U * is a ^-cover of X, 

i=l 

where \X\ denotes the diameter of X, that is, the maxi- 
mum possible distance between any elements of X: 

\X\ = sup{\x -y\:x,y e X}. 

Here, we say that a countable collection of sets Ui, with 
\Ui\ < 5, is a J-cover of X if X C Ug^E/*. 



2 



Notice that H also depends on a parameter 5 which 
expresses the scale under which the measure is taken. 
We can eliminate such dependence by applying a limit 
over <$, defining in this way the s-dimensional HausdorfT 
measure: 



H S (X) 



lim HUX). 



If we plot H S (X) as a function of s we observe a similar 
behavior in any fractal object analyzed. The value of H 
is oo for any s < D and it is for any s > D, where D 
is always a non- negative real value. D is the Hausdorff 
fractal dimension of X. In a more formal way we may 
write: 



When this dimension is estimated over a gray-level dig- 
ital image / : [M, N] — >> 5ft, an usual approach is to map 
it onto a three-dimensional surface S through: 

S = {i,j,I(i,j)\(i,j)€[l:M]x[l:N}}. (1) 

In this case, we construct a three-dimensional grid of 
3D cubes also with side-length S. The probability p m is 
therefore calculated as being the number of grid cubes 
containing m points of the surface divided by the maxi- 
mum number of points inside a grid cube. 



III. FRACTAL DESCRIPTORS 



D(X) = {s}|inf {s : H 3 (X) = 0} = sup {H S (X) = 00}. 

In most practical situations, Hausdorff dimension uses 
to be complicate or even impossible to calculate. Thus, 
assuming that any fractal object is intrinsically self- 
similar, we may derive a simplified version, also called 
similarity dimension or capacity dimension: 



D 



log(AQ 
log(r) ' 



where N is number of rules with linear length r used to 
cover the object. 

In practice, the above expression may be generalized 
by considering N as any kind of self-similarity measure 
and r as any scale parameter. This generalization gives 
rise to a lot of estimation methods for fractal dimension, 
with a broad application to the analysis of objects which 
are not real fractals (mathematically defined) but which 
present some degree of self-similarity in specific intervals. 
An example of such method is the probability dimension, 
used in this work and described in the following section. 



B. Probability Dimension 

The probability dimension, also known as information 
dimension is derived from the information function. This 
function can be defined in any situation where we have 
an object populating a physical space. We must divide 
this space into a grid of squares with side-length S and 
calculate the probability p m of m points of the object 
pertaining to some square of the grid. The probability 
function is given by: 
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where N is maximum possible number of points of the 
object inside a unique square. 

The dimension itself is given through: 



D 



lim 

5^0 



lnTVp 
TnT 



The main idea of fractal descriptors is to extract val- 
ues (descriptors) from the log — log relation common to 
most methods of estimating fractal dimension. Actually, 
any fractal dimension method derived from Hausdorff di- 
mension concept obey a power-law relation which may be 
explicit in the following: 



D 



log(aft) 

" log(e) ' 



where 9Jt is a measure depending on the dimension 
method and e is the scale under this measure is taken. 

Therefore, fractal descriptors are provided from the the 
function u: 

u : log(e) log(93T). 

In order to simplify the notation we name the indepen- 
dent variable as t. Thus, t = loge and our fractal de- 
scriptors function is denoted u(t). For the probability 
dimension used in this work, we have: 



u(t) 



\og{N v (5)) 
\og{5) ' 



The values of function u(t) may be used directly as de- 
scriptors of the analyzed image or may be post-processed 
by some kind of operation aiming at emphasizing some 
specifically aspects of that function. Here, we apply a 
multiscale transform to the function. In this way, we 
obtain a two-dimensional function U(b, a), in which the 
variable b is related to t and a is related to the scale on 
which the function is observed. A usual means of obtain 
U is through a derivative process: 

+00 

U(b, a) = ^ u(t — n)G(n, a), 

71— — OO 

where G is the well-known Gaussian function and a is 
the smoothing parameter: 



G(n, a) 



1 



:exp(-n 2 /2t). 



Figure 1. A diagram of probability dimension estimation. Above, the 2D version used for shapes analysis. Below, the 3D 
version used for gray- level images. 



Given the finite domain of function u, we also must re- 
strict the response of Gaussian filter to a finite interval 
[-M,+M]: 



+M 



U(b,a) = u(t — n)G(n,a), 

n= — M 

where M is a real value which should satisfy: 

roo roo 

2 G(t,a)dt = 2 G(v,a)dv<z, 

Jt=M Jv=M/y/a 

being e the tolerance error. Usually, we may have: 
M = Ccr + 1 = Cy/\a) + 1, 

where C is a real constant with values commonly varying 
between 3 and 6. 



IV. PROPOSED METHOD 

The idea here proposed is to obtain fractal descriptors 
from textures based on the probability fractal dimension. 
Thus, such descriptors are co mput ed from the curve u(t) : 
\og(Np(S)) in the Equation II B Therefore, we apply a 
multiscale transform to u. 

The multiscale process is achieved by deviating numer- 
ically and convoluting with a Gaussian filter, as described 
in the previous section: 




U(b,a) 



du 
~dt 



Figure 2. Multiscale transform. 



where G a is the Gaussian function descritized over the 
interval -10 < x < 10, x e Z. 

As multiscale transform maps a one-dimensional signal 
onto a two-dimensional function, it is a process which 
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generates intrinsic redundancies. We may find differ- 
ent approaches to eliminate such redundancies keeping 
only the relevant information^. Here, we adopt a simple 
method named fine-tuning smoothing in which U (6, a) 
is projected under a specific value ao of the Gaussian 
parameter. Here we tested values of a varying between 
0.1 and 5 and used that values which provided the best 
performance in the experiments. 

Finally, we selected a particular region from U(b,ao) 
to compose the descriptors. Empirically, we observed 
that the initial points in this curve are relevant to a good 
performance in our application. Then, we established 
a threshold t after which all points in the convolution 
curve are disregarded. Thus, the values in the curve 
U(i, ao), 1 < i < t are taken as the proposed descriptors. 




k 

Figure 3. Discrimination of texture image by the proposed 
descriptors. 



V. EXPERIMENTS 

In order to verify the efficiency of the proposed tech- 
nique, we applied the probability descriptors to the clas- 
sification of two benchmark datasets and compared to 
the performance of other well-known and state-of-the-art 
methods for texture analysis. 

The first classification was accomplished over the Bro- 
datz dataset, a classical set of natural gray level textures 
photographed and put together in a boolPH This dataset 
is composed by 111 classes with 10 textures with dimen- 
sion 200x200 in each class. 

The second data set is the Vistex, a set of color tex- 
tures extracted from natural scenes 24 . Here, we employ 
a version of the dataset in which we have 7 classes, each 
one with a variable number of images of 256x256 pixels 
and converted them to gray-level images. 

We compare probability descriptors to other 4 other 
techniques, that is, Gabor wavelets^, Co-occurrence 



matrix 26 , Gray Level Difference Method (GLDM]P, 
a multifractal approach described irP^ and Bouligand- 
Minkowski fractal descriptors^. 

Finally, we classify each descriptor by a hold-out pro- 
cess (half of data to train and the remaining to test) using 
K-Nearest Neighbor (KNN), with K = 1 and compare 
the results. 



VI. RESULTS 

The TableU shows the correctness rate in the classifica- 
tion of Brodatz dataset using the compared descriptors. 
The proposed method obtained the best result with a 
12% advantage. For this result we used a = 0.1 and 
threshold t = 17. A particular important aspect in Bro- 
datz data set is the reduced number of descriptors of the 
proposed approach. This point may be specially impor- 
tant in large data basis when the computational perfor- 
mance is more relevant. Furthermore, the small number 
of features avoid the curse of dimensionality, which prej- 
udices the reliability of the global result. 



Method 


Correctness Rate (%) 


Number of descriptors 


Gabor 


81.3 


20 


Co-occurrence 


53.5 


84 


GLDM 


52.2 


20 


Multifractal 


35.1 


101 


Bouligand- Minkowski 


47.6 


85 


Proposed method 


91.2 


17 



Table I. Correctness rate for Brodatz dataset. 



On the other hand, the Table [TT] shows the results for 
the Vistex textures. In this case, we obtained the best 
result by using a = 0.1 and t = 80. Again, the pro- 
posed approach provided the greater correctness with a 
2% advantage. Again, we have a good result in a data 
set which presents a lot of challenges once it is aimed at 
color analysis while the proposed approach is gray-level 
based. This aspect turns significant even a tiny classifi- 
cation enhancing. 



Method 


Correctness Rate (%) 


Number of descriptors 


Gabor 


85.7 


20 


Co-occurrence 


70.1 


120 


GLDM 


43.5 


20 


Multifractal 


20.1 


101 


Bouligand-Minkowski 


68.8 


85 


Proposed method 


87.7 


80 



Table II. Correctness rate for Vistex dataset. 



Finally, the Figures [4] and [5] show the confusion ma- 
trices of the methods with best performances. In such 
figures, we have the predicted classes in the rows and the 
actual ones in the columns. The number of classes in 
each configuration is given by the intensity of gray-level 
in each point (brighter points correspond to large number 
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of classes). In this kind of representation, a good descrip- 
tor must produce a matrix with a diagonal as brighter 
and continuous as possible and the minimum of brighter 
points outside the diagonal. In this sense, we see that, in 
Brodatz data, the probability descriptors presented ex- 
actly these characteristics, with almost no "hole" in the 
diagonal and with a lower density of brighter points out- 
side. 
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Figure 4. Confusion matrices in Brodatz dataset. a) Gabor. 
b) Co-occurrence, c) GLDM. d) Proposed method. 



In Vistex case, the matrices are not so distinguishable 
visually, even due to the similarity in the results. Thus, 
a perspective which we may use is analyze directly the 
number of corrected samples in each class. Based on this 
aspect, we observe in the legend bar that the proposed 
technique has its matrix normalized on a greater number 
of samples. This implies that although with a similar as- 
pect, our approach presented a higher number of samples 
classified correctly in each class. 

A general analysis of the results demonstrates that the 
proposed method overcame the compared ones in both 
data sets, using a reasonable amount of descriptors. Such 
results was expected from fractal theory once we have a 
lot of works in the literature showing the efficiency of 
fractal geometry in the analysis of natural textures. Ac- 
tually, fractal geometry presents a flexibility in the mod- 
eling of objects which cannot be well represented by Eu- 
clidean rules. The fractal dimension is a powerful metric 
for the complex patterns and spatial arrangement usu- 
ally found in the nature. Fractal descriptors enhance 
such ability providing a way of capturing multiscale vari- 
ations and nuances which could not be measured through 
conventional tools. More specifically, the probability de- 
scriptors here proposed aggregates a statistical approach 



to fractal analysis, composing a framework which sup- 
ports a precise and reliable discrimination technique, as 
confirmed in the above results. 



VII. CONCLUSION 

The present work proposed a novel method to extract 
descriptors based on fractal theory for texture analysis 
application. 

Here we obtained such descriptors by applying a mul- 
tiscale transform over the power law relation of fractal 
dimension estimated by the probability method. 

We tested the efficiency of the novel technique in the 
classification of well-known benchmark texture dataset 
and compared its performance to that of other classi- 
cal texture analysis methods. The results demonstrated 
that probability fractal descriptors are a powerful tool to 
model such textures. The method provides a rich way of 
representing even the most complex structures in texture 
images, being a reliable approach to solve a large class of 
problems involving the analysis of texture images. 
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Figure 5. Confusion matrices in Vistex dataset. a) Gabor. 
b) Co-occurrence, c) Bouligand-Minkowski. d) Proposed 
method. 
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